Quasiperiodic Bloch-like states in a surface wave experiment 
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Bloch-like surface waves associated with a quasiperiodic structure are observed in a classic wave 
propagation experiment which consists of pulse propagation with a shallow fluid covering a quasiperi- 
odically drilled bottom. We show that a transversal pulse propagates as a plane wave with quasiperi- 
odic modulation, displaying the characteristic undulatory propagation in this quasiperiodic system 
and reinforcing the idea that analogous concepts to Bloch functions can be applied to quasicrystals 
under certain circumstances. 
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The main difficulty towards the development of a sys- 
tematic analytic approach to the transport properties of 
quasiperiodic systems has been the absence of an anal- 
ogous Bloch theorem approach as used in the periodic 
case. In the first efforts to apply a modified version of 
the Bloch theorem, it was noticed that the dense spec- 
trum of quasiperiodic systems is dominated by only a 
few special reciprocal lattice points that may be taken 
to construct a quasi-Brillouin zone §. Thus, by consid- 
ering only the dominant Fourier components, the atomic 
distribution can be expanded in terms of a discrete aperi- 
odic lattice. Wave functions of the form \f/k = itk(r)e lk ' r 
will therefore solve the Schrodinger equation. In this case 
Mk is quasiperiodic and should formally be defined on a 
countable dense set of reciprocal lattice vectors. But, by 
the above considerations, this expansion is useful since 
the Fourier development of the modulation function Uk 
can be restricted to the few special reciprocal vectors that 
dominate the spectra. Thus, Bloch-like states could de- 
scribe the plane wave propagation in so schematized qua- 
sicrystals and free-electron-like bands are expected. Re- 
cently this idea was experimentally tested showing that 
analogous concepts to Bloch functions can be applied to 
quasicrystals 

The classic wave propagation in quasicrystalline sys- 
tems was addressed in a first seminal acoustic experi- 
ment of He and Maynard || by the feature that acous- 
tical waves are ideal tools to investigate formally similar 
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quantum propagation effects Q . On the other hand, ap- 
pearance of the quasicrystalline symmetry in fluids dy- 
namics was firstly predicted theoretically by Zaslavsky 
and co-workers Q and a simulation similar to the condi- 
tions of the present experiments was reported in Ref. || . 
Finally, compressible quasisymmetric flows were consid- 
ered in Ref. |t]], whereas a general outlook on order and 
disorder in fluid motion can be found in the experiments 
of Gollub §. 

In this letter we shall see that a discrete restricted 
spectral scenario can be displayed by means of impulsive 
waves in hydrodynamic quasicrystals, where we observe 
Bloch-like surface waves. The waves are generated at 
the frequencies corresponding to the Fourier components 
of the quasiperiodic structure at the dominant diffrac- 
tion spots. The observed Bloch-like waves are plane 
waves with quasiperiodic modulation generated when a 
pulse propagates transversally to the quasiperiodic struc- 
ture. Liquid surface waves shape a quasiperiodic grid 
that obeys the so-called Octonacci sequence, previously 
studied ^ |lO) but never observed in any experiment. 

The quasiperiodic structure involved in our experiment 
is the octagonal Ammann-Beenker tiling composed by 
squares and rhombuses jTTJ . Associated with this octag- 
onal tiling is a quasiperiodic sequence, named the Oc- 
tonacci sequence [|| , that can be generated starting from 
two steps L and S, which are related according to the ir- 
rational ratio L/S = 1 + V2, by iteration of substitution 
rules: L — > LSL and S — > L. 

The experiments are performed with surface waves 
generated on a shallow fluid that covers the quasiperi- 
odically drilled bottom of a transparent vessel. Such ex- 
periments are similar to others realized in vessels with 
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periodic bathymetry and described elsewhere |l2|, |13| but 
here not only a continuous wave excitation driven by a 
vertical monofrequency vibration but also a new experi- 
ment of transversal pulse propagation is performed. The 
bottom dimples are located at 121 vertices of the octag- 
onal tiling. The edge length I of the tiling is 8 mm with 
an error lower than 0.4%, the radius r of the cylindrical 
bottom wells is 1.75 mm with an error lower than 1% and 
their depth d is 2 mm. The depth of the liquid layer over 
the cylindrical wells is given by h 2 — hi + d, where hi is 
the depth of the thin liquid layer covering the bottom of 
the vessel among holes. 

Under conditions of continuous- wave excitation, an in- 
ertial hydrodynamical undulatory instability grows over 
the bottom wells when the system vibrates vertically at 
a frequency of 35 Hz. Such an instability becomes re- 
markable (Fig. |l|) due to the high density and the very 
low surface tension of the liquid |Q. Oscillating bulges 
over dimples are connected by surface waves with shorter 
wavelength that decorate the shallow liquid region among 
holes with hi being 0.4 mm. This is a physical scenario 
similar to that of the Kronig-Penney model but adapted 
here to a 2D quasiperiodic system. It should be remarked 
that Fig. |l| is the first available experimental example of 
a quasiperiodic pattern of waves not arising from a non- 
linear instability such as the Faraday instability. Stand- 
ing waves are coupled to only two rings of Fourier wave 
components of the quasiperiodic bottom structure, as 
shown in the inset of Fig. [l]. Their wavenumber ratio 
is 1 + \/2, as it can be measured in the diffraction pat- 
tern (see below). If the effective Fourier transform of a 
quasiperiodic structure is restricted to a discrete set of 
Fourier peaks as in this case, then Bloch-like modulation 
functions can be used to describe the wave propagation in 
such a simplified quasicrystal jll|. The pattern of Fig. |l| 
is due to the strong coupling between liquid surface waves 
and the bottom quasiperiodic topography and does not 
depend on the shape of the vessel boundary. Undistin- 
guishable patterns are generated with octagonal or cir- 
cular boundaries. Although the boundaries are reflect- 
ing vertical walls, the boundary symmetry matches the 
symmetry of this standing wave experiment allowing the 
mentioned boundary-independent strong wave coupling. 
Due to the accuracy on realizing the setup, localization 
phenomena do not appear in the described experiment. 
However, slight tiltings of the vessel generate wave do- 
mains [[12], pL3| ; furthermore, point and linear defects can 
be easily introduced in the system by dropping mercury 
on the bottom dimples to study new interesting wave 
localization phenomena |l6| ] . 

The propagation of a plane wave through the octago- 
nal quasiperiodic structure can be visualized by means of 
a experiment of wave pulse propagation. A coupling be- 
tween the vertical waves of the vessel-liquid system and 
the surface waves generated by a transversal pulse is ex- 
pected, modelling in this way the propagation of a plane 
wave through the quasiperiodic structure. This experi- 
ment was realized in a vessel with octagonal boundary. 



The octagon side L is 4 cm and it is perpendicular to the 
r — X direction of the well structure. The surface ratio 
/ between bottom holes and the whole octagon is about 
0.15 with hi being 0.5 mm. The system is excited near 
the octagonal boundary with a wave pulse parallelly to 
the liquid surface and perpendicularly to a side of the 
octagon. The signal is picked up by means of a Briiel 
& Kjaer 4344 accelerometer placed at the center of the 
vessel and it is processed by means of a digital acquisi- 
tion system. The impulsive signal and the corresponding 
Fourier transform are shown in Fig. ^ (inset and solid 
line, respectively). Three clear spectral peaks of the vi- 
brational vessel-liquid system appear at about 20, 30 and 
50 Hz. As we shall see, such resonances indicate the exis- 
tence of three narrow band gaps in the liquid surface wave 
propagation fl3[ ] , i. e. standing liquid waves are generated 
at approximately the above generated frequencies. 

At the start of each pulse the liquid feels the pertur- 
bation and a nice quasicrystalline surface wave pattern 
suddenly appears [Fig. ||(a)]. A transitory weak turbu- 
lence arises in the system after scarcely 0.04 s [Fig. ||(b)], 
whereas robust standing waves drawing clear quasiperi- 
odic grids can be observed between 0.08 and 0.24 s on 
the liquid surface [Fig. ||(c)]. Finally, quasiperiodic grid 
patterns decay until the arrival of the next pulse. As 
wave phase velocities are about 11 cm s _1 and the wave 
group velocity is nearly null near the gaps, times for an 
echo at the boundaries to come back are much longer 
than observation time. 

The robust quasiperiodically spaced standing waves 
shown in [Fig. ||(c)] are generated by discrete Bragg res- 
onances and thus can be considered quasiperiodic Bloch- 
like waves. To verify this, first note that the irrational 
ratio LS/L = \J~2 is apparent in our experiment [Fig. 
|(c)] . Now, using a crystallography-oriented computer 
program Jl7| , the Fourier transform of the pattern of 
Fig. ||(c) is calculated and shown in Fig. |] (top left). 
Such diffraction pattern matches with an adequate sub- 
set of the diffraction pattern of the direct product of both 
orthogonal Octonacci sequences calculated according to 
theoretical methods [[)[ [l^] as shown in Fig. |J (top right). 
The absence of some diffraction peaks indicates the di- 
rectional character of the impulsive action. The pulse 
runs along the T — X direction from the upper left to the 
bottom right corner in both patterns at the top of Fig. 
^. Along this direction, a intensity profile is taken in the 
experimental pattern and recovered the inverse Fourier 
transform of that unidimensional diffraction subset. The 
result is shown in Fig. ^ (bottom) which displays an 
Octonacci sequence, and it matches with that generated 
theoretically starting from the above mentioned substi- 
tution rules. Finally, the above described intensity pro- 
file along the T — X direction is scaled according to the 
wavenumber of the waves of the diffraction pattern shown 
in Fig. ||(c). Such scale is then changed according to the 
approximate dispersion relationship given p3[ by 

u? = gk { 1 + — k 2 ) tanh(£;/io), 
V P9 ) 
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where ho = ht(l — f)+hif, uj is the angular frequency, k 
is the wave number, g is the acceleration due to gravity, 
T is the liquid surface tension and p is the liquid density 
fl4| . In Fig. H (dashed lines) the gray scale intensity pro- 
file is plotted versus the frequency according to the above 
mentioned change of scale. The first maximum is scaled 
by the wavenumber k = 11 cm , that corresponds to 
the main ubiquitous wave appearing in the experimental 
pattern of Fig. ||(c) . As it can be seen, diffraction peaks 
which represent narrow band gaps closely match in fre- 
quency with those independently measured also in Fig. |^. 
Thus, in this restricted scenario, where the resonances of 
the vibrational coupling generates a discrete spectrum, 
the wave pattern observed in Fig. ||(c) corresponds to 
quasiperiodic Bloch-like states. 

If experiments of pulse propa gat ion are realized in ves- 
sels with periodic bathymetry |l2| , |i~3f no signal of tur- 
bulence appears. Thus, as remarked in a different con- 
text ||To| , 19[ , the quasiperiodicity of the hydrodynam- 
ical system could be the origin of the weak chaos ob- 
served in the described experiments just at the start of 
pulses, when amplitudes are higher and hence the non- 
linearity is stronger. Then, a rapidly increasing number 
of incommensurable Fourier harmonics can grow due to 
the finite frequency bandwidth of the pulse and the in- 
commensurate nature of the system. This gives rise to 
the pre-turbulent state of the surface waves. When the 



multiscattering becomes weaker, the Fourier mode cas- 
cade decays and the propagative wave exhibits a clean 
quasiperiodic grid pattern. 

Anyway, the quasiperiodic structure underlies in spite 
of the weak turbulence apparent in Fig. ||(b) , which must 
be looked at grazing incidence to recognize Octonacci 
quasiperiodic sequences. This is evident in Fig. ^, which 
is the inverse Fourier transform of Fig. ||(b). The direct 
product of two orthogonal Octonacci sequences is recov- 
ered there, showing a patch of the well known octagonal 
tiling filled with square and rhombic tiles ]il| . 

In conclusion, we have shown Bloch-like surface waves 
associated with a quasiperiodic structure in a classic 
wave propagation experiment. These waves draw clear 
quasiperiodic grids that obey the Octonacci sequence. 
Our results along with earlier ones H can be helpful to 
understand the characteristic undulatory propagation in 
quasiperiodic systems. 
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FIG. 1: Snapshot of the system vibrating vertically at a fre- 
quency of 35 Hz, under conditions of continuous wave excita- 
tion. The inset shows its Fourier spectrum with a well denned 
and discrete set of relevant components. 



FIG. 2: The solid line is the Fourier transform of the impulsive 
signal as picked up by an accelerometer at the center of the 
vessel. It shows clear resonances at approximately 20, 30 and 
50 Hz. Such resonances are specific of the system. The first 
peak at very low frequency, close to the origin, corresponds 
to the Fourier transform of the square pulse that excites the 
system. The signal in the time domain is presented in the 
inset. The dashed line is the gray scale of the subpattern along 
the r — X direction as represented in Fig. [j| versus frequency. 
The gray scale (between and 1) was rescaled to match the 
peak of the Fourier transform at 30 Hz. Three standing waves 
indicating narrow band gaps appear at approximately 20, 32 
and 50 Hz. 
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FIG. 3: Temporal sequence of patterns observed when the sys- 
tem is disturbed with a transverse pulse, (a) Quasicrystalline 
pattern observed at the start of each pulse, (b) A transitory 
weak turbulence is observed after 0.04 s. (c) Standing waves 
draw clear quasiperiodic grids between 0.08 and 0.24 s. This 
figure should be looked diagonally at grazing incidence. 





FIG. 4: Experimental (top left) and theoretical (top right) 
Fourier transform of the pattern shown in Fig. ^|(c). At the 
bottom, the inverse Fourier transform of a subpattern along 
the r — X is shown. The Octonacci sequence is clearly recov- 
ered. 



FIG. 5: Inverse Fourier transform of pattern shown in Fig. 
^|(b). The structure of the octagonal tiling underlying on 
the well quasiperiodic arrangement of the vessel bottom is 
recovered. 



